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I. INTRODUCTION 



In recent years modified theories of gravity of the f{R) type have been thoroughly stud- 
ied mainly because of the expected role that they could play in the understanding of the 
cosmic speedup. For different choices of the Lagrangian f{R), one can find cosmological 
self-accelerating solutions without the need for sources of dark energy [H, |2[. In this context, 
it was soon realized js], 0] that once an f{R) lagrangian is given, the field equations could 
be derived in two inequivalent ways depending on whether the connection is defined as 
the Levi-Civita connection of the metric (metric formalism) or whether it is seen as inde- 
pendent of the metric (Palatini formalism). The resulting equations are radically different 
from each other and only lead to the same dynamics in the particular case of f{R) = R — 2A. 



The main focus in the literature has been on f{R) theories in metric formalism, the 
usual variational formalism, which leads to higher order equations for the metric. Higher 
curvature terms of the f{R) type are well-known in the literature [sj and naturally arise in 
erturbative quantization schemes [6| and in low-energy effective actions of string theories 
which justifies the interest in this kind of theories. The Palatini formulation, however, 
was not known to make contact with more fundamental approaches and has only received 
a timid attention as compared with its nonidentical twin brother in metric formalism. 
This situation changed recently jsl, when it was shown that a certain Palatini Lagrangian 
could faithfully reproduce the effective dynamics of Loop Quantum Cosmology (LQC) |9[, a 
Hamiltonian based approach to quantum cosmology inspired by the non-perturbative quan- 
tization techniques of Loop Quantum Gravity [lOi] . Though the Hamiltonian description of 
metric f{R) theories has been discussed several times and from different points of view in 
recent literature [ll-13||, a discussion of the Palatini case is still missing. The main purpose 
of this paper, therefore, will be the study of the Hamiltonian formulation of Palatini f{R) 
theories. 



The dynamics of f{R) theories in metric formalism can be analyzed in several classically 
equivalent ways. Some authors treat them as higher derivative theories, whereas others in- 
terpret them as scalar-tensor theories. Using a covariant approach, one finds that the higher 
order derivatives that appear in the field equations always act on a particular combination 
of second-order derivatives of the metric, namely, the scalar curvature, R (or a function of 
it). Taking the trace of the metric field equations one then finds an independent equation 
which only involves (a function of) R and derivatives with respect to R of functions of R. 
This fact allows one to reinterpret the theory as a scalar-tensor theory instead of as an 
intrinsic fourth order theory, which simplifies the mathematics and makes the physics of the 
problem a bit more transparent. The scalar-tensor representation of the theory can be put 
into correspondence with a kind of Brans-Dicke like theories with characteristic parameter 
u = and a non-trivial potential, which is intimately related to the form of the function 
f{R): V{(j)) = Rdf/dR— f{R) and = df/dR. This technical advantage facilitated the 
analysis of cosmological models and of solar system tests. The results, however, are not very 



optimistic for those models in which f{R) is of the form R plus infrared corrections [14- 
Nonetheless, there is still hope that models with chameleon properties might be viable 
From a Hamiltonian point of view, it has recently been shown that the higher-derivative 
and the scalar-tensor interpretations are related by a canonical transformation 



|16|. 

a- 
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In the Palatini formulation of f{R) theories, metric and connection are seen as inde- 
pendent variables. Upon variation of the action, one finds that the connection satisfies a 
constraint equation, which can be easily solved in terms of the metric and the matter fields, 
whereas the field equations for the metric are second order, like in General Relativity (GR). 
The resulting theory has, surprisingly, the same number of degrees of freedom (dynamical 
fields) as GR. The role of the Lagrangian f{R) is thus to change the way matter generates 
curvature by adding new matter contributions, which stem from the connection constraint 
equation, on the right hand side of Einstein equations. In the case of infrared corrected 
models proposed to explain the cosmic speedup, such terms lead to incompatibilities with 
observations of many kinds [H. 13. [isl- 2l| . For ultraviolet corrected models, on the contrary, 
one finds the same agreement with observations as in GR for curvatures much smaller than 
the ultraviolet scale of the theory [12, 23 1. Furthermore, it has been shown that the 
equations of motion that follow from the effective Hamiltonian of Loop Quantum Cosmology 
for an isotropic cosmology can be derived from a covariant action of the f{R) type in Pala- 
tini formalism This cosmology turns out to be non-singular and replaces the disturbing 
Big Bang singularity of GR by a cosmic bounce, which occurs at a density of order the 
Planck density. This effect is seen as a success of the non-perturbative techniques of the 
loop quantization. The existence of a Lagrangian that reproduces the LQC dynamics is 
important not only because it establishes the covariance of that Hamiltonian theory but also 
because it shows that Palatini f{R) theories and its generalizations could be an appropriate 
framework to study aspects of the phenomenology of quantum gravity. More recently, it 
has been shown that bouncing cosmologies are quite generic in the Palatini f(R) framework 



that extended theories of the type f{R, Rfj_^R'^'^) also share that property j2J, |25|, and 



that fi^R) theories, where ^R is the scalar curvature of the Barbero-Immirzi connection 



with parameter /3, are classically equivalent to Palatini f{R) theories [26 



Palatini f{R) theories also admit a scalar-tensor representation a la Brans-Dicke 15l. |27 
In this case, however, the theory is characterized by the parameter u = —3/2 and a non- 
trivial potential, which has the same form as in the metric version of the theory . This 
particular value of u is very special because it leads to a non-dynamical scalar field, which 
justifies why the theory has the same number of dynamical fields as GR. In this sense, it is 
well-known that the existence of non-dynamical fields in a theory results in constraints in 
its corresponding Hamiltonian formulation. However, though the Hamiltonian description 
of generic Brans-Dicke theories is very well known, the particular case of Palatini f{R) 
{u = —3/2 and V{(f)) ^ 0) has not been studied yet. In addition, due to the interesting 
connections existing between Palatini j[R) and LQC, which is a purely Hamiltonian based 
theory, we find necessary a careful analysis of the Hamiltonian description of Palatini j[R) 
theories. A better understanding of the Hamiltonian theory will provide a new viewpoint 
for the comparison between the metric and Palatini versions of /(-R) theories and might also 
provide new insights on the possibility of finding effective descriptions for the dynamics of 
LQC and of more general Hamiltonian quantum gravity theories. 

Since the field equations of Palatini and metric /(-R) theories are (classically) equivalent 
to the Brans-Dicke cases w = —3/2 and = 0, respectively, in this paper we consider 
the Hamiltonian formulation of Brans-Dicke theories with a non-trivial potential rather 
than that of the original j{R) versions. This will allow for a more transparent comparison 
between metric and Palatini cases or, with more generality, between the case lo = —3/2 
and all the other Brans-Dicke cases u ^ —?>/2. We will show that the canonical momentum 
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associated to the Brans-Dicke scalar becomes degenerate in the case w = —3/2, whereas 
it remains an independent degree of freedom for any other value of the parameter w. The 
w = —3/2 case, therefore, must be handled with constrained system jlsj. We will 

see that its constraint and evolution equations share many similarities with the w ^ —3/2 
case though it also has some important differences. We will use these similarities and dif- 
ferences to comment on the well-jposedness of the Cauchy problem for w = —3/2, which has 



been criticized in the literature [2l, l29|] (see also [30|, |31[, and [32| for a different opinion) 



Being mathematically rigorous, we would like to remark that for the original Palatini 
f{R) one should consider the 3 + 1 decomposition of the metric and of the connection inde- 
pendently, which complicates the analysis and the comparison with the metric version of the 
theory. However, all the applications of Palatini theories found in the recent literature (in- 
cluding its applications to LQC and nonsingular cosmologies) are restricted to the properties 
of their equations of motion and, therefore, can be translated to the Brans-Dicke case with- 
out further discussion. For these reasons, at this stage we find an unnecessary complication 
the explicit 3-1-1 decomposition of the original Palatini f{R) theory. Such decomposition, 
however, seems unavoidable in more general Palatini theories whose Lagrangian contains 
other scalars besides R (such as R^^R^^ , R^lvapR^'"^^) for which an on-shell scalar-tensor 
representation is neither known nor likely to exist. 

The content of the paper is organized as follows. In section //, we present the covariant 
formulation of f{R) theories and describe their relation with Brans-Dicke theories. We then 
use the scalar-tensor representation to work out the 3+1 decomposition of the action previous 
to the construction of the Hamiltonian. The next step is to compute the Hamiltonian for 
general Brans-Dicke theories with w ^ —3/2 and compare our result with the literature. 
The case w = —3/2 is worked out and compared with the general case w ^ —3/2 in section 
V. We derive, compare, and discuss the constraint and evolution equations in Sec. VI, 
where we also comment on the Cauchy problem. We conclude with a brief summary and 
conclusions. 



II. COVARIANT FORMULATION 

Let us begin by defining the action of Palatini f{R) theories 

S[g,r,^^] = d^x^fXR) + Sm[g,u.i^m] (2.1) 

Here f{R) is a function oiR = g^'R^u{T), with R^,{T) given by R^,{T) = -d;r\^ + dxTl,+ 
^^p^ux ~ ^up^^^J.\ where F^^, is the connection. The matter action Sm depends on the matter 
fields ipmi the metric g^^, which defines the line element ds^ = g^^dx^dx^ , and its first 
derivatives (Christoffel symbols). By construction, we assume that the matter action does 
not depend on the connection F^^,. Varying (12. ip with respect to the metric g^^p we obtain 

fRR,u{V) - \f{R)g,. = (2.2) 

where fji = df /dR. From this equation we see that the scalar R is algebraically related to 
the energy-momentum tensor via the trace equation 

RfR - 2/ = K^T, (2.3) 
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The solution to this algebraic equation will be denoted hj R = R(T). The variation of (12. ip 
with respect to T'^^ must vanish independently of (12 ■2p and gives 



(2.4) 



where fR = fii{R[T]) must be seen as a function of the matter terms by virtue of (12. 3p . This 
equation is equivalent to Va [\/— = and leads to 

r;. = -J- (5mV + dutp, - 9,v) (2.5) 
where t^j, = fng^u- Expressing the connection in terms of fn and g^i,, (12. 2p becomes 

JR ^JR 

3 . . . 1 . ,n1 

+ 



2(/k) 



+ -J- [V,.V,/r - g^uUfn] 

jR 



(2.6) 



One should note that all the R, /, and terms on the right hand side of this system of 
equations are functions of the trace T of the matter energy-momentum tensor. In vacuum, 
T = 0, those equations boil down to 



R^lu{g) - ^g^,uR{g) = -A^jfg^u 



(2.7) 



where Ag// = (-Ro/i?o ~/o)/2/i?o is evaluated at the constant value Rq = R(T = 0) and plays 
the role of an effective cosmological constant. 



A. Equivalence with Brans-Dicke theories 



The field equations (12. 6 p derived above can be rewritten in a more compact and illumi- 
nating form by introducing the following definitions 

= fR (2.8) 
V^(0) ^ /?(0)/k - fiRm (2.9) 

where represents a scalar field and V{(j)) its potential. Note that in eq. (l2.9p we have 
assumed invertible the relation between R and /ij(-R) to obtain R{fR) = R{(p)- The equations 
of motion for the metric can then be expressed as follows 

1 1 

Rf^Ag) - ^9f.uR{g) = -jT^^ - —g^uV^ct)) + 



df,(j)d^(j) - ^g^u{d(j)f 



+ 



+ 7 - (?^.n0] (2.10) 
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where u takes the value u = —3/2. It is straightforward to verify that f{R) theories in 
metric formahsm lead to a similar set of equations but with the replacement u = (see for 
instance [Isj for details). In that case, the trace equation is SD//? + Rfn — 2/ = k^T. 

The equation of motion for the scalar field in both metric and Palatini formalisms is 
provided by the trace of the original field equations and can be expressed as 



dV 

(3 + 2w)n0 + 2V(0) -0— 



k'^T 



(2.11) 



Remarkably, the field equations f l2.10p and (12. lip can be derived from the following action 
S[g^,,<P,^^] = d^x^[ct>R{g)- (2.12) 



UJ 



- -{d,<pdy) - 1^(0) 



which represents a Brans-Dicke like scalar-tensor theory (in the usual metric variational 
formalism). Note that in the original Brans-Dicke theory the potential term was absent. 



V{(j)) = [33|. The absence of a potential term leads to inconsistencies when oo = —3/2, since 
then (12. lip becomes ill defined unless all matter sources satisfy T = 0. For non-trivial V{(j)), 
that equation makes perfect sense, since then can be expressed as an algebraic function 
of the trace, = 0(T), which becomes constant in vacuum (and wherever T =constant). 



III. HAMILTONIAN FORMULATION 

Having obtained the (on-shell) equivalence between f{R) theories and Brans-Dicke theo- 
ries, we use the scalar-tensor representation to work out the 3 + 1 decomposition of the action 
needed to construct the corresponding Hamiltonian. With regard to the field equations, this 
procedure leads to the same classical results as with the original f{R) theory. We choose to 
work with the equivalent scalar-tensor action because the f{R) representation in the Palatini 
formalism would require considerable extra effort for the implementation of the 3 + 1 slicing 
of the space-time and the analysis of the constraints jl^] due to the a priori independence of 
the connection. In the scalar-tensor form, the independent Palatini connection has already 
been factored out and appears in the form of a scalar field coupled to the curvature and 
characterized hj u = —3/2 and a non-trivial potential V{(j)). We will see that this specific 
value of the parameter u is interesting in its own right, though we prefer to maintain its 
relation with f{R) theories to emphasize the motivations that led us to study this particular 
case. 

From now on we use lower-case latin letters to represent space-time indices. 



A. 3+1 decomposition 

The 3 + 1 decomposition of the action fl2.12p proceeds in the usual way j35[. We consider 
a foliation of the spacetime manifold Ai into hypersurfaces St of simultaneity characterized 
by a function T{x) =constant, a normalized timelike co- vector oc daT normal to this 
hypersurface, and a shift vector N"- orthogonal to n'^ = g''''nb. This allows us to construct a 
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time flow vector = Nn°- + A^", where is known as lapse, and decompose the metric in 
the form gab = hah + snarly. The parameter s = ±1 = na'nf' tells us whether the signature is 
Euclidean, s = 1, or Lorentzian s = — 1. Elementary, though lengthy, manipulations allow 
us to express the Lagrangian density of fl2.12p as follows 

\/h 



N^h'''Da(pDb(p + s 



N^Da^ 



-2K (j) - N''Da<p) - NV (0) 

where Kab = hahfV is the extrinsic curvature. Da* 
scalar of the 3— metric hob, and we have used the following relations 

^(4) _ ^(3) _ ^ \^KabK''^ - {Ka^'f + 2^ J^] 

NrVc(l) = sh'"^D^(l)DdN + K (^(p- N^Da 



V\9\ 

The canonical variables of the theory are {gab, 4>) = {N, N"", hab, 
are deflned by the following expressions 

SS 

0, 



jjab 



ss 

^ = , n„ 

6N 
6S 



(3.1) 

t"9a0, (^^i? is the Ricci 

(3.2) 
(3.3) 

(3.4) 
(3.5) 

). The canonical momenta 

(3.6) 
(3.7) 



2^2 



SS 



J^ab 

h {K"^ - Kh"^) - — 



2k^ 



2K+^{^-N^D^, 



Like in GR, we immediately see that the momenta conjugated to and N°- are con- 
strained to vanish. On the other hand, from the combination of Uh = hab^""^ and tc^, we 
flnd that 



3 + 2u; 
N 



(3.9) 



is also constrained to vanish when u = —3/2, which is the case associated with Palatini 
f{R). At this point, it is useful to rewrite the Lagrangian density C using the deflnition for 
to eliminate the explicit dependence on Kab from it. The result is 

S (2/«2^2 



Nlo 



-D^<i)D''(t) + 2D^<pD''N - NV{<p) 
+ (1 + s)2K U - N'D,(j)) + (3 - 2sw) 



(3.10) 



N^Dr 



21 



2Ncf> 
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The first term in the last line of this expression vanishes in the Lorentzian case, s = — 1, 
and the second term also vanishes when s = —1 if u = —3/2. When (3 + 2w) ^ 0, the 

combination (j) — N'^Dc(j)j in the last line can be expressed in terms of the momenta using 

(13. 9p . To proceed with the construction of the Hamiltonian one must have in mind the above 
constraints and apply Dirac's algorithm for constrained Hamiltonian systems. From now on 
we take s = —1. 



B. General case uj / —3/2 



In the general case cu ^ —3/2, we have the same (primary) constraints as in GR, namely, 
Ctv = IlNit,x) = and Ca = Ila{t,x) = 0. The Hamiltonian is constructed by introducing 
Lagrange multiplier fields X]sf{t,x), \a{t,x) for the primary constraints and performing the 
Legendre transform as usual with respect to the remaining velocities. The result is 



H = j d^x [XnCn + X^Ca + N^-Ha + NUn] 



where 




H. 



{2k 



2\2 



hcj) 



n2 

Wu,, - ^ 1 + 



-De0^>+V(0) + 



{2k 



2\2 



2/i0(3 + 2uj) 



(n, 



"Ha 



(3.11) 

(3.12) 
(3.13) 

(3.14) 



These expressions reproduce previous results found in the literature for Brans-Dicke theories 
with V^(0) = [37]. For the dynamics to be consistent, the constraints must be preserved 
under evolution, which requires that = {H, Cjy} = and Ca = {H, Ca} = 0, where the 
poisson bracket at time t is defined as 



{A{x),B{x')} = / £a 



6A{x) 5B{x') 5B{x') 5A{x) 
(a) 5Qi{a) ~ 5Ii\a) 5Qi{a) 



(3.15) 



where H' and Qi generically represent the canonical variables. In particular, this definition 
leads to {W\t,x),h,d{t,x'] = 5^,5^)5(2^ 

X — x'). By direct evaluation, one finds that Cn = 

—5H/5N = —T-Ln and Ca = —SH/6N°- = —T-La- We thus see that on consistency grounds 
we must impose the secondary constraints 71^ = and T-La = 0, which implies that the 
Hamiltonian H is constrained to vanish, like in GR. If matter is present, one must add the 
corresponding pieces 5Hmatt/5N and 6Hmatt/SN"- to these constraints, which leads to 



(3) 



H^ 



+ 2h"'D,Dd(j) + Vi(P) + 



2 



(3.16) 



2DdK + ^<^Da(p + 



20(3 + 2uj) 

1 ST-Lmatt 



+ 



15% 



matt 



a SN 



a (5iV« 







(3.17) 
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where we have defined a = h^l'^ / {2^?) and used the tilde to denote the tensorial quantities 
= '^4>/ot and 11"* = Il^^/a. The evolution equations can be written as follows (this requires 
some lengthy calculations which we omit here) 



N 



3 + 2w 



dV 



+ -D^D'^ct) - — 



+ 2wD^ [ ^D^(p 1 - 2AN + N^Dc^^ 



N{Ilh-(Pfc^)inh-2(P^^) 
202(3 + 2w) 

2N 

2 



+ 



hab = 2Di^aNb) + —\Ilab-^-^hab\ + 



0(3 + 2w) 



l-ab 



tl"^ = -N 



+ - f w^ul - — n™"n. 



2 



+ D''D\N(P) - h''^A{N(f)) - 2D^ND^(j) + h^^D^ND 

_ ^/.'^f _ 1 ^Hmatt _ 5iV(nfe-07f^) ^,, 
2 a 6hab 20(3 + 2u;) 



(3.18) 



(3.19) 
(3.20) 



(3.21) 



where we have denoted h"-^DaDhf = A/. The equations for vf^ and 11"* have been obtained 
from those corresponding to tt^ and 11"* by using the generic relation vr = vr/a — 7fg"*gaf,/2, 
where tt = {H, vr}. Using the tensorial quantities vr instead of the densities vr, we remove the 
term a from the constraint and evolution equations everywhere except in the contributions 
coming from the matter Hamiltonian. 



C. Palatini f{R) case: uj = -3/2 

As we saw above, the case u = —3/2 has an additional constraint not present in the 
general case u ^ —3/2. Equation (13.91) puts forward the fact that the momentum 'k^ of the 
Brans-Dicke scalar can be expressed as a linear combination of the momenta associated to 
the 3— metric hab- Taking the Lagrangian f l3.10p particularized to w = —3/2 and following 
Dirac's algorithm, the Hamiltonian of this case becomes 



(3.22) 
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where 



H 



N 



Vh 

2^ 



D^<PD^(P + V{<p) 



(3.23) 
(3.24) 

(3.25) 



and the other constraints are the same as in the general case w ^ —3/2. The evolution 
equations are now given by 



2AN + 2wD^ [ ^D^A + N^Dan^ - - 



hah = '^D(aNb) + ^ ( Ilafc - ^'^h ) + >^<phab 



1 

i 

2 



2 a 5hab 2 



(3.26) 
(3.27) 



(3.28) 



(3.29) 



Consistency of the evolution implies the following secondary constraints for = and 

Ca = 



+ 2h"'D,Dd(j) + V{(j)) + 



1 5H 



matt 



a SN 







-2D^K + n,Da<P + l'-^ = 0. 



(3.30) 
(3.31) 



Using the evolution equations and the constraint ( I3.30p . one can verify that the evolution 
of Cw, leads to 



— {H, C^} 



(3.32) 



O AUrrA\ I ATA^^ ^ ^'Hmatt , SHmatt 
2al\V[(p) + aN(p—— — hab- 



2 5N 



5h 



ab 
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Since must vanish, we must impose the secondary constraint 

Using the definitions T„fe = --^^jp^^ and ^f"^ = h"^ - ^ - N") [t^ - N^) to derive the 
relations 

5C _ 5£5g^ _ -2 , _ 5C_ 

5N ~ 5g^^ 5N ~ N^^ " ^ bg^^ ^^'^^^ 
h-^^ = r^i^ (3.35) 

§}^ab ^gab V J 

T — ^ j —N 6Cmatt ^ab ^^matt \ /„ q(-\ 

and the fact that = and ^fp^ = -^^1^, one can verify that fl3.33p yields 

0^ - 2\/(0) = K^T . (3.37) 

This equation reproduces the relation (12. lip when w = —3/2 and establishes an algebraic 
relation between the trace of the energy-momentum tensor of matter and the scalar field 

= 0(r). 



IV. ON THE EVOLUTION AND CONSTRAINT EQUATIONS 

Besides the number of constraints, the only difference between the Hamiltonians (I3.1ip 
and f l3.22p is the last term appearing in the definition of in (I3.14p . which is proportional 
to {Ufi — 07r^)^. That term stems from the last factor of (I3.10p . which obviously vanishes 
when u = —3/2 and, therefore, cannot appear in (I3.25p . For w = —3/2 one must introduce 
the constraint to account for the degeneracy imposed by the vanishing of (Ilh — (pn^)- 
This small difference in the Hamiltonians means that the constraint and evolution equations 
will share many terms in common. We have arranged the formulas in the text in a way 
that allows for a clear comparison between the w = —3/2 and the general w ^ —3/2 case. 
We thus see that the constraints Cn = 0, Ca = 0, and Ha = are the same in both cases, 

whereas Hn = and ^Hn = differ only slightly. Eq. (l3.16l) contains the term ^20(3+2^) ' 
which is absent in (I3.30p . As a result, the relation between the evolution equations can be 
summarized as follows 

N /- ^ \ 

(f)(n,) = + X<t>4>- J^^^y^h- (f>T^cf>j (4.1) 

<p ^ 2 * ^ 202 3 + 2w ^ ^ 

^1:^ = h[-b'^ - ^^Kb + (n. - 0^<^) Kb (4.3) 
n?:, = n- J, + -A,n- - (4.4) 



12 



These relations make it clear that there is no single (on-shell) choice of the function A,^ that 
allows to recover the w ^ —3/2 equations from the w = —3/2 case, unless one makes the 
obvious substitution — )■ 2h(p{3+2uj) ~ ^'^'t>) directly in the Hamiltonian f l3.22p . 

Now that we have compared the structure of the constraint and evolution equations, 
we will comment on their physical meaning and implications. In the case w ^ —3/2, the 
dynamical variables are (0, vr,^, /lah, 11"^) (plus the {qi,p^) of the matter), while the lapse 
and shift, (A^, Na), although necessary to determine the evolution, are undetermined by the 
equations. Note also that they do not enter into the constraint equations, for a change on 
the lapse-shift pair leaves the fields at the initial surface unchanged, which is a manifesta- 
tion of the freedom to specify the coordinate system. In the w = —3/2 case, the dynamical 
variables are just {hab, H"'') (plus the {qi,p^) of the matter), because the evolution equations 
for (0, TT^), as we saw above, can be combined to establish the secondary constraint fl3.37p . 
In this sense, it is worth noting that the constraint f l3.16p involves up to second-order spatial 
derivatives of hab and (see the terms ^^^R and DaDb(j)), but only first order time deriva- 
tives of them (contained in the momenta 11"'' and tt^). However, though the w = —3/2 
constraint fl3.30p contains spatial derivatives of = 0(T) up to second order, it does not 
contain any time derivative of 0(T) because the corresponding momentum tt^ is absent in 
that equation. Something analogous occurs in the vector constraint f l3.3ip . where we can 
use the replacement tt^ = n;i/0 to show that no extra time derivatives of the matter appear 
in the constraints. This is a very important aspect, because it means that the highest order 
time derivative of the matter fields appearing in f l3.30p and fl3.3ip is the same as in GR 
and coincides with the highest order present in the energy-momentum tensor of the matter. 
The evolution equations also have this property. A glance at 03.26113.2^ puts forward that 
the evolution equations for 0, n^, hab, and 11"^ do not contain the momentum vr^. Therefore, 
though one can find up to second-order spatial derivatives of 0(T), and hence of T, there is 
no trace of extra time derivatives acting on the matter fields. 

It is worth noting that the second-order spatial derivatives of 0(T) require an extra de- 
gree of smoothness in the matter profiles, an aspect that is not necessary in GR. This extra 
degree of differentiability is a natural requirement if we attend to the f{R) formulation of 
the w = —3/2 theory. Since the affine connection is compatible with a metric tab which is 
conformally related with the space-time metric Qab, the smoothness and differentiability of 
the conformal geometry is guaranteed if the conformal factor is differentiable up to second 
order (to yield a smooth field strength, Riemann tensor, of the affine connection). Since 
the conformal factor is the function /ij(T) = 0(T), the differentiability condition on the 
geometry is transferred to the higher-order (spatial) differentiability of the matter fields 
living in that geometry. This should not come as a big surprise because modified theories 
of gravity usually lead to higher-order equations for the metric, which demand a higher 
degree of differentiability of the metric field. This, in particular, occurs with f{R) theories 
in metric formalism. In the Palatini case, we do not find higher-order derivatives of the 
metric because the modified dynamics is due to the new role played by the matter fields. 
However, the existence of a conformal geometry intimately related with the matter fields 
ends up imposing the extra (spatial) differentiability conditions on the matter fields. 
Now, the existence of up to second-order spatial derivatives but not of time derivatives of 
{4>,TT^) when w = —3/2 is an unsolved question that requires further stud}0. In this sense. 



In fact, since the field (j) is not dynamical, the time evolution of the pair ((^, tt,/,) must be completely 
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we believe that a clearer understanding of the algebra of constraints in the case w = —3/2 



could shed some useful hght into this problem [38 



V. ON THE CAUCHY PROBLEM 

It is well-known that if in GR one specifies initial values for A^, A^"^, hab and 11'^'' which 
are consistent with the constraint equations, the evolution equations uniquely determine 
hab and 11"*, while N and A^** remain undetermined, which expresses the existing gauge 
freedom of the theory. This guarantees that the intrinsic (coordinate-independen t) g eome- 



try of space-time is determined uniquely by an initial choice of hab and 11'^'' [35|, |39| . The 
same is true for the scalar-tensor theories considered here, thus implying that the initial 
value problem is well- formulated for all w. For the w = —3/2 case, the only difference with 
respect to GR is that one must specify an initial value for and take also into account its 
corresponding constraint equation to consistently establish the initial data. 
Using different variables and representations of the evolution and constraint equations, one 
can also proof the well-posedness of the initial value problem of GR and of generic Brans- 
Dicke theories with w ^ —3/2 in both Einstein and Jordan frames ji^]. One can also make 
special choices for the lapse-shift pair and manipulate the corresponding 3 + 1 equations 
of GR to show that the conjugate variables hat and 11"* do satisfy a hyperbohc evolution 
system j4l|. Although we will not try here to explicitly demonstrate the well-posedness 
of the initial value problem for w = —3/2, we will exploit the resemblance between our 
constraint and evolution equations (I3.30l) . (l3.3ip . (l3.37p and (I3.28p - (l3.29p with those of GR 
to argue that the Cauchy problem is likely to be well-posed also for the Brans-Dicke case 
w = -3/2. 

Note first that in vacuum, T^^, = or Hmatt = 0, the constraint (13.370 implies that is a 
constant, 0o, which turns the constraints (13.300 and (13.310 into 

- <PoR^'^ + ^ (fl'^'flab -^]+ V(0o) = (5.1) 



^0 



2 



- 2D,Ul = . (5.2) 

With a simple constant rescaling of the metric, these constraints are the same as those of 
GR with a cosmological constant. Setting for consistency the Lagrange multiplier = 0, 
the evolution equations for hab and 11"* also recover the same form as those of GR with 
a cosmological constant. We can thus conclude that the Cauchy problem in vacuum is 
well-posed. 

When matter is present, one should add to our system of equations those corresponding 
to the matter fields, which we assume standard. The strategy now would be to interpret 
the 0-dependent terms, which are functions of the trace T, as part of a new (or modified) 
matter Hamiltonian. This way, the constraint and evolution equations maintain a structure 

provided by the time evolution of {qi,p''), as is manifest from the secondary constraint p.37p . One 
could thus expect that the possible higher-order time derivatives of the matter fields induced by time 
derivatives of (j) or tt^ in the field equations could be rewritten in terms of first-order time derivatives 
using the evolution equations for the matter variables (qi,p^). However, as we are remarking here, no such 
higher-order time derivatives arise and, therefore, we do not need to worry about this problem. 
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that closely resembles that of GR except by some non-constant factors 4>{T) that multiply 
or divide objects like ^^^R and n^^flaft. If the matter fields satisfy the spatial differentia- 
bility requirements imposed by the constraint equations, the absence of higher-order time 
derivatives of the matter fields in the constraint and evolution equations suggests that the 
time evolution will be as well-posed as in GR. This, in fact, has been explicitly shown for 
a perfect fluid using the Einstein frame representation of the evolution equations [32[. Ob- 
viously, since in general the well-posedness of the GR equations depends on the particular 
matter sources considered, the modification of the source terms induced by the existence of 
0(T)-dependent terms requires a model by model analysis. Therefore, though one cannot 
conclude that the Cauchy problem is well-posed for an arbitrary f{R) Palatini Lagrangian, 
we find no reasons to suspect that it is ill-posed in general. 



Before closing this section, we would like to mention that the Cauchy problem has been 
used in recent literature to criticize the viability of all Palatini f{R) theories. In fact, in [i^ 
it was claimed that the disappearance of the d'Alambertian 00 from (12. lip for the value 
w = —3/2 implies that the non-dynamical field can be arbitrarily assigned on a region or 
on the entire spacetime, provided its gradient satisfies a degenerate equation [Eq. (4.5) in 
that paper], which reduces to a constraint. This fact, it was stated, would make impossible 
to eliminate the term 00 from the evolution equations unless 00 = 0. This was interpreted 
as a no-go theorem for Palatini f{R) gravity, which would have an ill-formulated Cauchy 
problem even in vacuum. This interpretation is conceptually wrong because the scalar field 
in the w = —3/2 is just a given algebraic function of the trace T and, therefore, is clearly 
specified by the matter content. Moreover, one should note that Eq.(4.5) of 29| is not cor- 
rect. That equation should recover the well-known relation 21^ — 01^' = k^T that establishes 
the algebraic relation between and T [our secondary constraint (I3.37P ]. Using Eqs. (3.4) 
and (3.5) of [29], it is easy to check that the associated Eq. (3.10) does recover our equation 
fl3.37p in the Brans-Dicke case w = —3/2 (even though this is not the result obtained in 
^29i] ) . This indicates that the first claims against the well-posedness of the Cauchy problem 
for Palatini f{R) theories stemmed from the analysis of erroneous equations. 
The strong conclusions of [2^ were a bit relaxed in [2| (see in this sense j30|, [3l[ ) , where it 
was admitted that the Cauchy problem should be well-posed in vacuum and with radiation 
fields (for which T = and =constant). In fact, in [2^ it was correctly noticed that in 
the w = —3/2 case the field could be algebraically solved in terms of T (though their Eq. 
(219) is the same as Eq.(4.5) of (29|). It was then argued that the existence of terms of 
the form 00 (T), which imply contributions of the form DT, would cause problems for the 
Cauchy problem. Though such terms and the possible existence of higher-order derivatives 
of the matter fields are certainly a reason for concern, it was prematurely concluded that 
the Cauchy problem for Palatini f{R) theories was likely to be neither well-formulated nor 
well-posed unless the trace T were constant. These conclusions contrast with our findings, 
which show that the evolution equations do not introduce higher-order time derivatives of 
the matter fields, which guarantees that the initial value problem is as well formulated as 



in GR 35, 39 
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VI. SUMMARY AND CONCLUSIONS 



In this work we have carried out the Hamiltonian formulation of Brans-Dicke theories 
with a non-trivial potential paying special attention to the case w = —3/2, whose field 
equations are equivalent to those of an f{R) theory in the Palatini formalism. We have 
found that the scalar field in the case w = —3/2 presents a degenerate momentum, which 
is proportional to a linear combination of the momenta of the induced metric hab- This 
degeneracy requires, following Dirac's algorithm for constrained systems, the introduction 
of a new constraint in the Hamiltonian. Consistency of the evolution of that constraint 
leads to a secondary constraint which establishes an algebraic relation between the scalar 
field and the trace of the energy-momentum tensor of the matter. We have written the 
constraint and evolution equations in a way that allows for a quick comparison between the 
general case w ^ —3/2 and w = —3/2. 

We pointed out that the resulting constraint and evolution equations of the case w = —3/2 
do not contain any higher-order time derivative of the matter fields, and only spatial deriva- 
tives of the scalar 0(T) appear up to second-order. This implies that the spatial profiles 
of the matter sources must satisfy stronger differentiability requirements than in GR. We 
have interpreted this property as a natural requirement due to the existence of a confor- 
mal geometry directly related with the matter fields. By comparing the constraint and 
evolution equations of our theory with those of GR we have shown that the initial value 
Cauchy problem is well-formulated in general (because the intrinsic geometry of space-time 
is determined uniquely by an initial choice of hab, H"'', plus the corresponding position and 
momenta of the matter on the initial Cauchy surface) and have found reasons to believe that 
it is likely to be also well-posed, which contrasts with other opinions found in the literature. 
This comparison is possible because the equations do not contain time derivatives of either 
the scalar 0(T) nor of its momentum conjugate ir^. The reason for the absence of time 
derivatives of these objects is an unsolved question that may be related with particular 
properties of the enlarged algebra of constraints of the case w = —3/2, which will be studied 
elsewhere SSl . 
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